Incoherent bremsstrahlung in flat and bent crystals 
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Abstract. Incoherent bremsstrahlung by high-energy particles in crystal is due to the thermal 
spread of atoms in relation to their equilibrium positions in the lattice. The simulation procedure 
developed earlier for the incoherent radiation is applied to the case of the electrons and positrons 
motion in the sinusoidally bent crystal. The results of simulation are in agreement with the 
data of recent experiments carried out at the Mainz Microtron MAMI. The possibility of use of 
the sinusoidally bent crystals as undulators is discussed. 



1. Introduction 

The bremsstrahlung cross section for relativistic electrons in a crystal is split into the sum of 
the coherent part (due to the spatial periodicity of the atoms' arrangement in the crystal) and 
the incoherent one (due to the thermal motion of atoms in the crystal) [U [3]. Although 
the spectrum of incoherent radiation in crystal is similar to one in amorphous medium, the 
incoherent radiation intensity can demonstrate substantial dependence on the crystal orientation 
due to the electrons' flux redistribution in the crystal (channeling etc.). The simulation based 
on the semiclassical description of the radiation process [H [6] confirms that viewpoint. The 
results of simulation are in a good agreement with the corresponding early [7] and recent [8] 
experimental data. 

Here we present the results of simulation using the improved procedure taking into account 
the crystal deformations. The simulation was carried out under the conditions of the recent 
experiment performed at the Mainz Microtron MAMI [9] to explore the radiation emission from 
periodically bent crystal. The possibility of application of such crystals as undulators is discussed 
during last years [91 [TOl E] . 

For the reader convenience, in the next section we outline some theoretical ideas of our 
approach. 

2. Computation method 

Let us consider the high energy electron incidence on the atomic string in the crystal. The 
two-dimensional multiple scattering angle # is equal to the sum of individual scattering angles 
on the atoms: 

= J>(Pn). (1) 

n 

where p n is the impact parameter of the collision with the n-th atom of the crystal. The mean 
square of the multiple scattering angle (averaged over the thermal vibrations of atoms) can be 



expressed as the sum of two blocks of terms: 
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The first block describes the coherent scattering which can be interpreted as a motion in the 
uniform string potential j2[ [31 [12]. The second one describes the incoherent scattering of the 
electron on the thermal vibrations of the lattice atoms. 

The bremsstrahlung spectrum of the electrons passing through the crystal also can be 
expressed as the sum of the coherent and incoherent contributions [TJ [2j [3] . For the electrons of 
the energy e ~ 1 GeV the main contribution of the coherent effect is made to the soft range of 
the spectrum (the photon energy ftw is less or of the order of dozens of MeVs) . In the medium 
and hard ranges of the spectrum the incoherent radiation is predominant. 

The spectral density of the incoherent radiation from the individual electron moving on the 
given trajectory is described by the formula [6] 



d£ 
du J 



where m and e are the electron's mass and charge, c is the velocity of light. It is convenient to 
compare the incoherent radiation intensity of the uniform electron beam in the crystal to the 
corresponding intensity in the amorphous medium (described by Bethe-Heitler formula). The 
ratio of these two values is equal to [5[ [6] 
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where N is the total number of the electron's collisions with atoms under its motion through 
the crystal, R is Thomas-Fermi radius of the atom, integration over d 2 pQ means the integration 
over all possible points of incidence of the electron on the crystal in the limits of one elementary 
cell. 

The values of the function F(p) = (#(p) 2 ) — (#(p)) 2 are determined by linear interpolation 
of the values pre-calculated on the regular grid of impact parameters. The impact parameters 
p n are found using the electron's trajectory obtained by numerical integration of the equation 
of motion in the field of the set of parallel uniform strings. The influence of the incoherent 
scattering by the thermal vibrations of atoms on the electron's trajectory can be taken into 
account by adding to each component of the electron's transverse velocity the random value with 



the dispersion c^J |(#(p) 2 ) — (#(p)) 2 | /2 after each collision. For the further computational 
details see [5j [6] . 



3. Results of simulation 

The origin of the orientation dependence of the incoherent radiation intensity is illustrated on 
the figure 1. When the electron is incident to the atomic plane of the crystal under angle 6 less 
than some critical angle 9 C , it can be captured by the attractive potential of the plane. The 
finite motion in that potential is called as the planar channeling (see, e.g., [2"1 131 [12]). 

Under planar channeling the electrons collide with atoms under small impact parameters 
more frequently than in amorphous medium, that leads to the increase of the incoherent 
bremsstrahlung intensity; for above-barrier motion the situation is opposite. The account of the 
incoherent scattering of the particles on the thermal vibrations of the atoms leads to dechanneling 
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Figure 1. (a) Typical trajectories of the electrons and positrons under planar channeling. 
Pluses mark the positions of atomic strings (perpendicular to the figure plane) forming the 
atomic planes of the crystal. The horizontal scale of the figure is highly compressed, (b) The 
same for above-barrier motion, (c) Simulated incoherent bremsstrahlung intensity (in ratio to 
the Bethe-Heitler intensity in amorphous medium, see (dJ) from 1 GeV electrons vs incidence 
angle 9 to (110) plane of 30 fim thick Si crystal [I]. Dotted line corresponds to the trajectories 
simulated neglecting thermal vibrations of atoms, (d) The same for positrons 

and, hence, to the smoothing of the describer orientation dependence (compare solid and dotted 
lines on the figure 1 (c) and (d)). 

In the present article the simulation was carried out under the conditions of the experiment 
[9], where the radiation from e = 855 MeV electrons under their incidence onto the silicon 
crystal with sinusoidally bent (110) planes had been studied. The yield of photons with the 
energy hw = e/2 (for which the incoherent mechanism of bremsstrahlung is predominant) had 
been registered. 

The results of simulation (figure 2) demonstrate the qualitative agreement with the 
experimental data [9J. We can see characteristic structures similar to one on the figure 1 (c), 
generated by different crystallographic planes with the common [110] axis. 

The decrease of the radiation intensity in comparison to the reference flat crystal permits 
to estimate the reduction of the dechanneling length due to the crystal bending. We can see 
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Figure 2. The incoherent bremsstrahlung intensity (in ratio to Bethe-Heitler intensity in 
amorphous medium) from 855 MeV electrons in flat (upper plot) and sinusoidally bent (lower 
plot) silicon crystals under scanning of the goniometric angle <j) like in the experiment [9] 



that the bending of the crystallographic planes increases the dechanneling rate so highly that 
the scattering on the thermal vibrations of atoms already have no substantial influence on the 
incoherent radiation intensity (compare solid and dashed curves on the lower panel of the figure 
2). 

The main cause of the rapid dechanneling in the bent crystal lies in the arising of the 
centrifugal addition to the planar potential [12] : 

x 

U e ff = U(x) — e— under \x\ ^> Rb (5) 
Rb 

(where Rb is the bending radius) and, as a consequence, to the reduction of the potential 
barriers between which the channeling could take the place. For the crystal [9] the bending 
radius Rb > 6.2 • 10 -3 m. Assuming the planar potential has the shape of quadratic parabola 
with the depth Uo ~ —23.5 eV for (110) plane, the effective potential © in the domain of 
maximal curvature of the crystal would have the shape presented by the curve 2 on the figure 3. 
We see substantial decrease of the potential wells depth comparing to the reference flat crystal 
(curve 1). So, only a small part of incident particles could be captured into the channeling 
regime. 



Figure 3. Potential energy 
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Figure 4. The same as on the figure 2 for positrons. 



4. Conclusion 

The incoherent bremsstrahlung yield exceedance (deficiency) from the electron (positron) beam 
in the oriented crystal in comparison to the value predicted by Bethe-Heitler formula for the 
amorphous medium can be used as an indicator of the relative amount of the particles moving 
in the channeling regime. 

The results of simulation demonstrate that the planar channeling effect leads to the 
orientation dependence of the incoherent bremsstrahlung intensity not only in flat, but also in 
bent crystals. However, in the bent crystal the electrons rapidly leave the planar channels. The 
figure 4 demonstrates that for the positron beam instead of the electron one the dechanneling 
rate would be almost the same. The effective capture of particles into planar channels (and 
so the use of sinusoidally bent crystals as undulators) is possible under changing the bending 
parameters to ensure smaller curvature than in the crystal [S]. 
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